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I. INTRODUCTION 


Tht; research program is organized Into three different ■subjects: 
(1) the Jupiter system, (2) the Saturn system, and (3) comets. In the 
Jupiter system, the extended atmospheres of lo and their relations 0 ;.d 
interactions with the planetary magnetosphere are the topics of inter- 
est. In the Saturn system, the hydrogen torus of Titan, the atmospheres 
associated with the planetary rings, and the extended atmospheres > f the 
E-Ring satellites are the topics of concern. For comets, the dust and 
gas clouds in the outer coma are the area of focus. 

Research efforts for the extended atmospheres of lo are divided 
into the study of the lo sodium cloud, the lo potassium cloud, and other 
gas clouds of lo. The effort discussed here represents a logical con- 
tinuation, extension and expansion of our prior programs supported by 
NASA. A condensed summary of our prior and prenent programs is given 
in Table 1. As is shown in Table 1, the study of the lo sodium cloud 
is further divided into four separate topics organized around different 
types of observational data that are available. Efforts in the first 
year of this piogram have been priamrily focus*;- upon model improvement 
or development for the topics 1 (i) , l(iii) and 2 in Table 1, and upon 
exploratory investigation for the topic 1 (iv) . The topic l(ii) has 
remained inactive, waiting for the first year improvements in the model 
to be completed and the topic 3(ii) will be activated when relevant new 
information becomes available. 

The research effort for the Saturn system also represents a logi- 
cal continuation and an extension of our prior programs supported by 
NASA. A brief summary of the content of the prior and present programs 
is given in Table 2. As can be seen in Table 2, the Titan hydrogen 
torus has been a long-standing topic, whereas interest in the extended 
ring and E-ring satellite atmospheres has been initiated more recently 
as a result of the Pioneer 11, Voyager 1 and Voyager 2 encounters with 
Saturn. For this reason, effort in the first year of the program has 
placed primary emphasis upon topic 2, where some data are now available, 
and minor effort on topics 1 and 3. These latter two topics will re- 
ceive increased emphasis in the work to be performed in the second and 
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and Application to 
Observations 



third years (see Appendix A for a brief post-Voyager update on the Titan 
hydrogen torus) . 

For comets, the research program here represents a continuing of 
effort to refine our basic exoaphe'iic model for comets used for inter- 
pretive analysis of both dust and neutral gaseious components in the 
coma. The prior and present programs for this development are also sxjm- 
marized in Table 2. No effort to continue this development was sched- 
uled during this first program year in order that more emphasis could be 
placed upon the model developments and improvements scheduled under 
topics l(i) and l(iii) of Table 1. The necessary refinements and the 
subsequent application of the comet model are to be suitably scheduled 
during the next two years. 

In Sections II, III and IV, the progress and significant accom- 
plishments of the first year of the program are sun arized. References 
are given in Section V and Section VI, and are followed in Section VII 
by the Appendices containing supporting information. 
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11. PROGRESS IN MODELING THE 10 SODIUM CLOUD 


2. 1 Model Development 

The primary focus of our effort during the first year for the lo 
sodium cloud has been the generalization of the cloud model to include 
the time and spatial dependent sink produced by the oscillation of the 
lo plasma torus about the satellite plane. Electron Impact ionization 
of sodium atoms by torus electrons provides the onl sink mechanism con- 
sidered to date, but elastic collision of sodlun cloud atoms and plasma 
tc. us ions ma; also possibly contribute (Brovm, Pilcher and Strobel, 

1981) . The Importance of the elastic collisions will be evaluated in 
the second year effort. 

The oscillation of the plasma torus about the satellite plane 
introduces several different and distinct time dependences into the 
sodium cloud for an earth observer. Relative to the satellite, the 
period of oscillation of the plasma torus is about 13 hours and this 
will produce u time dependent sink for the cloud with observed intensity 
changes on the order of a few hours. Such changes have been observed and 
attributed by Trafton (1977, 1980) to the oscillating torus. Apparent 
time dependent changes in the two-dimensional morphology of the cloud 
also occur on a longer time scale for an earth observer because of the 
periodic motion of lo about its circular orbit ('^2.5 hour) and the dif- 
ferent geometrical projections of the cloud Intensity on the sky plane. 
Because of the non-synchronous values of the oscillatory and orbital 
periods, the sodium cloud on the sky plane will actually appear to have 
a longer effective period of about 25 lo rotations or A4.25 days. Real 
time dependent changes having this same 44.25 day period will also be 
introduced because of the east-west cloud distribution asymmetry which 
is driven by solar radiation pressure (Smyth, 1979, 1982). These real 
time changes occur because the envelope of the sodlun cloud changes 
shape as a fuiction of the satellite phase angle and therefore has a 
different spatial overlap with the oscillating plasma torus for each 
different angular location of the satellite on its orbit about Jupiter. 
This different spatial overlap is illustrated for two diametrically oppo- 
site satellite pliase angles in Figure 1 and is described in more detail 
in Appendix B. 
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The generalization of the sodluzs cloud otodel to include Che 08 Cll~ 
lacing plasma torus requires processing the time history of many atom 
orbits as they travel through the oscillating torus. This time proces- 
sing, although relatively straight forward, involves careful alteration 
and considerable expansion of the irevious existing computational code 
with special emphasis placed upon execution efficiency* Steady and sub- 
stantial progress has been made and the first phase, that of model devel- 
opment to li^cliide Che oscillating torus excluding the effects of solar 
radiation pressure, has been successfully ccx<tpleted. The second phase, 
that of adding the effects of solar radiation pressure to the general- 
ized code structure created in Che first phase, was Just completed at 
the end of this first program year. The improved model is capable cf 
calculating the two-dimensi'inal density and intensity distributions of 
the sodium cloud on the sky plane [topic l(i) of Table 1] as well as 
the shape of the spectral line profile seen through a prescribed view- 
ing aperture [topic l(ii) of Table 1]. The model is also capable of 
calculating the instantaneous and thirteen-hour cumulative spatial dis- 
tribution of sodium ions deposited in the planetary magnetosphere. 

Very preliminary results for the model developed mder phase one 
have been obtained during code refinements and code testing. A ”one- 
orbit-run" option was mostly used in these tests instead of the full set 
of orbits necessary to provide the complete spatial statistics for the 
cloud. This allowed the code to be thoroughly tested while at the same 
time minimizing the CPU computing expenses incurred. Fluctuations by a 
factor of two or more in the cloud density north and south of To were 
obvious in these test runs and were correlated with the spatial location 
of the oscillating torus. The average lifetime of sodiua in the cloud 
increased significantly because of the oscillating torus anu provided an 
increased abundance of sodiimi in the cloud for a given excape flux at lo 
when compared to the non-oscillating computation. Sodium ions that were 
created were preferentially deposited near the intersection of the satel- 
lite orbit plane aud the centrifugal plane of the plasma torus. In the 
beginning of the second program year, preliminary results will be cal- 
culated for the full model completed under pliase two of the development. 


12 



2.2 lo Plasma Torus 


Efforts to Improve the description of the lo plasma torus in the 
lo sodium cloud model have been undertaken simultaneously with the model 
development summarized in the previous sub-section. The lifetime of 
sodiun atoms in the plasma torus due to electron impact ionization is 
presented in Figure 2. The calculation is based upon the best avail- 
able infcrtoatlon from analysis of Voyager 1 data. Electron densities 
derived fr<^m plasma data by Bagenal and Sullivan (1981) and electron 
temperatures deduced both from UVS data by Shemansky (1980) and from 
in situ plasma measurements by Scudder, Sittler and Bridge (1981) were 
used. 

In Figure 2 , the lifetime of sodium atoms is seen to be as small 
as about one hour near lo and to increase at lo's location during the 
7 degree oscillation of the plasma torus about the satellite plane to 
a value of about 4 hours. Along the centrifugal eqiiator, the lifetime 
of sodium is extremely asymmetric about lo's orbital location. This 
asymmetric character will allow atoms which travel inside of lo's orbit 
to have a larger survival lifetime than atoms that travel outside of 
lo's orbit. This difference in survival lifetime may provide a basis 
for explaning the previously assumed asymmetric inner-hemisphere source 
of sodiun used to match observational data in model calculations where 
the plasma torus was not included. Future model calculations to be 
performed in the second year will explore this possibility more fully. 

Efforts to extend the accuracy of the plasma torus description in 
the model, especially for radial distances larger than 8 Jupiter radii, 
have also been pursued during the first year. Voyager 1 in situ plasma 
data describing the complete radial region from 4.9 to 21 Jupiter radii 
have been obtained from Belcher (1981) and are under study for this rea- 
son. A suitable description for the plasma density in this extended 
radial region has been adapted from thi.s study and will be incorporated 
in the model in the second program year. Ground-basid observations 
(Pilcher et al. , 1981; Trauger, 1981) showing changes in the plasma 
properties with the magnetic coordinate will be incorporated in 

the future when a definite fuictional dependence can be established. 
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Sodiuoi Electron Impact Ionization in the lo Plasma Torus. 
The lifetime calculation is based upon the plasma data 
referenced in the text. The two-dimensional common tem- 
perature model of Bagenal and Sullivan (1981) was selected 
to describe the electron density. 





2.3 Directional Features 


The peculiar directional features, discovered by Pilcher (1980) 
and also observed by Goldberg (1981), represent one of the newest phen- 
omena of the sodium cloud to be documented. From the data of P.''.lcher 
(1981), these directional features are seen to be somewhat transient In 
nature and are usually observed outside of lo's orbit as faint exten- 
sions to the bright sodium cloud. These features many times appear to 
be inclined about 20 degrees with respect to the satellite plane and 
sometimes appear above and sometimes appear below the satellite plane. 

Preliminary modeling analysis of Pilcher's data for these 20 
degree inclined features, undertaken in the first year, suggests that 
high speed sodlisn (6-10 km/sec) is emitted near lo with a velocity vec- 
tor most likely in the forward direction of. motion of the satellite and 
has a similar inclination angle with respect to the satellite orbit 
plane. Analysis of more data will be required to determine if this 
high-speed emission is localized on the satellite surface, or if the 
emission results from the magnetospheric wind that moves past lo. In 
the latter case, such a high-speed stream of sodium might be driven by 
elastic collisions of sodium atoms with plasma torus ions (Brown, Pilcher 
and Strobel, 1981) and provide a source for the remote sodiun cloud ob- 
served outside of the lo orbit by Brown and Schneider Q981). 

The possible magnetospheric wind emission process is illustrated 
in Figure 3 and contrasted with an inner hemisphere emission process. 

The inner hemisphere emission process was required for much lower vel- 
ocities ('''2.6 km sec ^) in earlier models of Smyth and McElroy (1978), 
which predated the discovery of the lo plasma torus, to obtain reason- 
able agreement withthe two-dimensional morphology of the bright sodium 
cl )ud observed by Murcray and Goody (1978). The actual initial velocity 
distribution at lo is likely broad (0-18 km sec ^) with a peak value 
probably at the lower emission velocities. One method being studied to 
help determine the actual relative weights in tha initial velocity dis- 
tribution is that of using the sodium model including the ionization 
effects of the plasma torus to study data (Pilcher, 1981) for the time 
evolution of the directional features. Such a study may also be import- 
ant in evaluating if there is a variable flux of sodium escaping from 
lo's exosphere which is modulated through the magnetospheric wind inter- 
action by the oscillatory motion of the plasma torus about the satellite. 
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Figure 3. Emission of Sodium from lo. Two rather different schemes 
by which sodium may be emitted from lo are Tustrated. 



111. PROGRESS IN MODELING THE 10 POTASSIUM CLOUD 


J.l Introduction 

' “ O 

The lo potassium cloud has been observed to radiate In the 7665 A 

O 

and 7699 A lines by Trafton (1975a). The radiation is thought to arise 
from solar resonance scattering. Observations of Trafton (1977» 1981a) 
have shown that the lo potassium cloud exhibits spatial and temporal 
variations similar to those observed for the sodium cloud. These obser- 
vations indicate that the potassium cloud isidergoes periodic fluctuations 
in response to solar radiation pressure (producing the east-west asymmet- 
ries) and in response to the ionizing Influence of Jupiter's plasma torus 
(producing the north-south asymmetries). In addition, temporary direc- 
tional features of potassium from lo have also been observed. Modeling 
improvements for the sodium cloud model discussed in the preceeding sub- 
section may therefore be applied directly to our future development of 
an lo potassium cloud model. 

Efforts in the first program year tiave been scheduled to obtain 
the necessary information for potassiizn required for converting, in the 
second program year, the lo sodiisn cloud model to the lo potassiiin cloud 
model. This information consists of the electron impact ionization life- 
time of potassium in the lo plasma torus, the solar spectrum for each 
potassium line, and the acceleration experienced by potassium atoms as 
they undergo solar resonance scattering. All this information has been 
acquired in the first program year and is simuoarized below. 

i. 2 lo Elasma Torus 

The electron impact ionization lifetime of potassiixn in the lo 
plasma torus has been calculated and is presented in Figure 4. The 
calculation is based upon the same Voyager 1 data used in determining 
the sodium lifetime in Figure 3. On the centrifugal equator, the potas- 
sium lifetime is slightly smaller near lo's orbital location and has a 
less steep gradient with increasing radial displacement from Jupiter 
than does the sodiun lifetime. The ionizing impact of the plasma torus 
on the potassium cloud may, therefore, be expected to be slightly more 
severe than for the sodium cloud in our future modeling. 
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3. 3 Solar Spectrum for the Potasslua Lines 

Both of the potar^ium lines appear in solar absorption features, 

O 

and the line at 7665 A is almost completely obscured to a ground-based 
observer because of an optically thick oxygen absorption feature present 
in the earth's atmosphere. To properly model the intensity morphology 

O 

of the cloud or the line profile shape in the 7665 A emission, the exact 
shape of the solar absorption feature above the earth's atmosphere must 
be knovm. The computation of the unobstructed solar absorption feature 
and the earth's atmospheric transmission function has been provided in 
this program year by Kurucz (1982) with whom a collaborative effort for 
this purpose has been established. The results of his computation are 
presented and discussed below. 

The calculated solar spectrum above the earth's atmosphere, the 
calculated transmission function of the earth's atmosphere at the Kitt 
Peak Observatory and the calculated and measured solar spectra at the 
location of the Kitt Peak Observatory are given in Figure 5 and Figure 

O O 

6 for the potassium lines at 7699 A and 7665 A respectively. The cal- 
culated and measured spectra are in excellent agreement and differ at 

O 

most by 1%- For the potassium line at 7699 A in Figure 5, the atmo- 
spheric transmission is essentially 100% except for the two very small 
absorption features in its far wings. The depth of the potassium fca- 

O 

ture at 7699 A is about 19% of the continuum value at line center and 
the wings of the feature above the earth's atmosphere are essentially 

O 

symmetric. For the potassium line at 7665 A in Figure 6, the earth's 
atmospheric transmission function is zero or small over much of the 
solar absorption feature. Above the earth's atmosphere, the potassium 
feature has a depth of about 14% at line center and the wings of the 
feature are asymmetric about the line center. This asymmetry results 

O O 

because of weak CN solar absorption features at 7664.316 A, 7664.404 A, 

o o 

7664.533 A and 7664.550 A and stronger Fe solar absorption features at 
7664.162 A, 7664.305 A and 7664.411 A. 

The physical implications of the solar spectra and the earth's 
atmospheric transmission functions presented above may be translated 
into interesting effects in the observed and modeled spatial and spec- 
tral characteristics of the lo potassium cloud. Two different types 
of effects should be distinguished. First, there are brightness varia- 
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cions introduced because the intensity contribution available for solar 
resonance scattering by a potassium cloud atom depends upon the shape 
of the solar spectrun and the amount of Doppler shift along this spec- 
trum produced by the relative radial motion between the atom and the 
sun. Second, there will be additional brightness impediments for a 
ground-based observer which are introduced by the shape of the earth's 
atmospheric transmission feature and the amount of Doppler shift along 
this transmission spectrun produced by the Earth-Jupiter velocity and 
the motion of the cloud about Jupiter. These two effects are discussed 
below. 

The radial velocity between a potassium cloud atom and the sun 
varies from nearly zero to a value as large as 20 km sec ^ to 35 km 
sec ^ and will provide a Doppler shift from line center as large as 
500 ma to 875 ma. Doppler shifts of only about 250 ma are required for 
the potassium atom wavelength to be shifted from its line center to its 
wing. Larger Doppler shifts can only provide a significant change in 

O 

the intensity of the 7665 A spectrum, and then only for Doppler shifts 
larger than 400 ma to shorter wavelengths (when the potassium cloud is 
observed near western elongation) because of the presence of the strong 
Fe absorption feature in the solar spectrun (see Figure 6). The radial 
velocity bet\^een the potassium cloud and the earth may vary from nearly 
zero to a value as large as 40 km sec ^ to 65 km sec ^ and will provide 

O O 

Doppler shifts as large as 1.0 A to 1.625 A relative to the wavelength 

O 

scale of the earth's atmospheric transmission function. For the 7665 A 

O O 

line, Doppler shifts between 0.75 A and 1.25 A to a larger wavelength 
will be located in the next adjacent oxygen absorption feature (see 
Figure 6) so that care must be taken in comparison of observed and cal- 
culated spatial intensities or high-resolution line profile shapes to 
properly account for the atmospheric transmission function. 

3.4 Solar Resonance Scattering 

O 

From Figure 5 and Figure 6, the potassiun lines at 7699 A and 

o 

7665 A are seen to be located in solar absorption features with depths, 
respectively, of about 19Z and 14Z of the continuum level and with widths 

O 

of about ±250 mA aroiind the line centers. This width is sufficiently 
narrow that potassium atoms with about 10 km sec ^ of radial velocity 


22 



respect to the sun have emission that are Doppler shifted completely 
out of the absorption feature. In fact* for a potassiim atom moving 
about Jupiter with lo's circular speed of 17.32 km sec continuum 
level sunlight is available for resonance scattering over more than 
half of its orbital path. This is in contrast to a sodltm atom in 
Io*s extended atmosphere which has an absorption depth of about 5X of 
the continuum level and an absorption width sufficiently broad that 
the full contiuuim level is note reached even for a Doppler shift cor- 
responding to a radial velocity of 17.32 km sec ^ with respect to the 
sun. 

This difference in the shapes of the potassium and sodium solar 
absorption features will be reflected directly in a different functional 
dependence, with repsect to Che satellite phase angle, of the cloud in- 
tensity and Che perturbing effects of solar radiation pressure on the 
cloud. In contrast to this different functional dependence, the maxlmun 
magnitudes of the acceleration produced by the solar radiation pressure 

will in fact be comparai e for both potassiim and sodium atoms, having a 

•> -1 

value of about 1.25 cm s and 1.36 cm sec re..pectively for a 17.32 
km sec ^ Doppler shift. These two accelerations correspond respectively 
to 1.7% and 1.9% of the gravitational acceleration of Jupiter at lo's 
orbit. In modeling, therefore, the lo potassiim cloud should exhibit 
east-west asymmetries of similar magnitude to those discussed earlier 
for the sodiian cloud by Smyth (1979, 1982). 



IV. PROGRESS IN MODELING SATURN'S RING ATMOSPHERE 


4 . i Introduction 

A hydrogen atmosphere associated with Saturn's rings has been 
reported from rocket measurements (Welser, Vltz and Moos, 1977; Welser 
and Moos, 1978), from the earth-orbltlng Copernicus (Barker et al., 

1980) and lUE (Clarke et al., 1981) satellites and from the Pioneer 11 
(Judge, Wu and Carlson, 1980) and Voyager 1 ind 2 (Broadfoot et al., 
1981; Holdberg, 1982) spacecrafts. Analysis of the rocket data by 
Carlson (1980) suggested that photodlssoclatlon of surface water mole- 
cules on the Icy ring material might be the source of hydrogen for the 
atmosphere. A preliminary survey of Information describing charged- 
particle sputtering and photon-sputtering of H 2 O Ice and the chemical 
produc'^s so formed was therefore undertaken to further explore this 
idea . 

in this survey, virtually no information has yet been uncovered 
to characterize the photo-sputtering process. Measurements for electron 
and ion sputtering on icy surfaces have, however, been pioneered by L. J 
Lanzerotti, W. L. Brown and colleagues at Bell Laboratories In collabor- 
ation with R. E. Johnson at the University of Virginia (see Brown et al. 
1980 ). Applying these measurements to the Saturn sytem, Cheng, Lanzer- 
otti and Pirronello (1982) concluded that the relevant ion phase space 
densities show that the dominant contribution to charged-particle sput- 
tering arises in the vicinity of Saturn's E-rlng and the moons Dione and 
Tethys. The phase space density above the A and li rings is too small to 
provide a significant source of hydrogen for the ring atmosphere. The 
authors also noted that the photo-sputtering rate assuned by Carlson 
( 1980 ) for the nearly maxlmun tilt angle of the rings with respect to 
the sun, appropriate for the rocket measurements of Welser, Vltz and 
Moos ( 1977 ), should decline significantly when the rings are nearly 
edge-on to the sun, a geometry appropriate to the Voyager observations. 
The results of the rocket measurement and Voyager UVS measurements are, 
however, comparable, which led the authors to question the validity of 
the photo-sputtering mechanisms. 
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Preliminary efforts in the first year hiave been initiated to 
address such apparent discrepancies. A study has therefore been mder- 
taken to investigate the sources, sinks and orbital dynamics of H atoms 
for the hydrogen ring atmosphere as well as the absorption, re-emlssion 
and chemical nature of H 2 O icy surfaces. This is viewed as a prelimin*- 
ary and necessary step before developing a serious nwaerical model for 
the ring atmosphere. As an initial departure point, effort has been 
mainly restricted to studying (1) the initial velocity distribution of 
hydrogen expected from photo-sputtering and (2) the fraction of the lib- 
erated hydrogen that is spatially restricted to the near vicinity of the 
rings. 

4 . 2 Preliminary Study 

In the absence of available information for photo-sputtering, a 
preliminary study was begun to look at the upper limiting energetic 
case for photo-sputtering of H 2 O ice: that of photo-absorption of H 2 O 

vapor by the solar UV radiation. Thi should be the upper limit case 
since all binding energy of H 2 O molecules with the ice surface and all 
inelastic absorption of energy by the ice surface are neglected. Results 
of our survey of this energetically limiting upper case are summarized 
in Table 3. The results are divided into reactions for four wavelength 
intervals and the product velocities for H atoms and OH molecules are 
indicated- The fraction of the total H 2 O photoabscrbed for each reac- 
tion is also given for both quiet and active sun conditions. For quiet 
sun conditions, 90. 7X ot the reactions produce H atoms dire tly, 3.4% 
of the reactions produce H 2 molecules and 0 atoms directly, and 5.9% of 
the reactions produce various ion products. For the 90.7% of the reac- 
tions producing H atoms, 67% have speeds peaked at about 20 km sec ^ 
from reaction la, 16.7% have speeds of 30 km sec ^ from reaction 2a, 4% 
have speeds of about 5 km sec ^ from reactions 2c, 2d, and 3c, 2.8% have 
speeds varying between 25 tc 35 km sec ‘ from reaction 3a, and 0.2% have 
speeds between 0 and 17 km sec ^ from reaction 3d. From this it can be 
concluded that 95.4% of the H atoms directly liberated have velocities 
between about 20 and 35 km sec ^ (2,09 to 6.4 ev), that 4.4% have vel- 
ocities of about 5 km sec ^ (0.13 ev) or less and that 0.2% have veloc- 
ities between 0 and 17 km sec ^ (0 to 1.5 ev). 
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rable 3. Photoab&orption of Solar UV Radiation by H^O Vapor 
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9S-A H O + hv — ioniiition products (ll 0 , etc.] — -- 6. 9 .059 .069 



The A, B, C and D rings of Saturn are located radially between 

abou: 1.0 and 1.3 planetary radii (R ). For a nominal ring radius of 

® -1 
1.3 R , the circular orbital speed of a ring particle Is 20.5 km sec 

and It has a period of 2.73x 10 sec. For an H atom emitted from the 
ring particle to escepc from Saturn's gravitational field, It must have 
a minimum viloclty in the forward direction of motion of 29.0 km sec ^ 
or a velocity vector relative to the ring particle of 8.5 km sec ^ 

(i.e., a kinetic energy of 0.38 ev) . For an H atom emitted frcMn the 
ring particle to collide with Saturn, It must have a mlnlmun velocity 
in the backward direction of motion of only 2.17 km sec ^ (I.e., a 
kinetic energy of 0.025 ev) relative to the ring particle. It is then 
clear from the results of Table 3 above that most H atoms produced from 
water vapor will be lost from the rings by collision with Saturn and 
g'-avitationa] escape from Saturn. In fact, a simple calculation reveals 
that the fraction of H-atoms that have Initial orbits within the radial 
interval between one and four Saturn radii is only 7Z. Subsequent col- 
lision with the ring particles of this selected 7% will cause even these 
K atoms to be rapidly lost through collisions with Satuin. The residence 
time of an H atom in the torus would under these c ire urns tar.ces be too 
short to allow a substantial H cloud to be formed for the rings. 

If an effective binding snergy is assumed lor the water ice sur- 
face to describe the work function of the surface and all Inelastic ab- 
sorption of energy by the ice surface, the results described above for 
the pure water vapor case can be generalized and significantly improved. 
._viiy state balances for the nvxnber of H atoms in the rings can in fact 

.established which allow residence times of sufficient lengths to pro- 

33 

. u the required 8 x 10 H atoms in the ring atmosphere measured by 
Weiser et al. (1977) and Weiser and Moos (1978). Even though this ap- 
peals o be a promising solution for formation of the H atom ring atmo- 
spher' , c her more difficult questions such as how to limit the abun- 
dance of the atomic oxygen ring atmosphere (produced by the photoabsorp- 
tion reactions of Table 3) below the upper limit of 2.6 Rayleighs set 

o 

for the 1304 A lines by the UV measurements of Voyager 1 (Broadfoot et 
al., 1981) remain to be answered. This and other related issues will 
continue to be pursued In our further studies during the second program 
year . 
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Introduction 


Slide 1 


I will report on modeling of observat ms for UV and 
optical emissions from extended neutral gas atmospheres of 
the outer satellites. 

Three objectives of this modeling are to better under- 
stand the local satellite atmosphere, to better understand 
the interaction of the satellite and its atmosphere with the 
magnetosphere, and to better assess the importance of the 
extended satellite atmosphere as an ion-source for the mag- 
netosphere . 

The presentation today will focus only upon the last of 
these objectives by discussing modeling results for the ion 
creation rate and the initial spatial distribution of the 
satellite ion-sources produced by lo's atomic oxygen cloud 
and by Titan's hydrogen torus. 

IP's Oxygen Cloud 

The atomic oxygen clnud of lo was discovered by Brown 
in 1980 by earth-based observations of its emission in the 
6300A oxygen lines. 

Our model calculation for the lo oxygen cloud viewed 
at the mid-point of Brown's measurements is shown in the 
first slide. 

Brightness contours for the 6300A emission are indicated 

in Rayleigh and assume a satellite escape flux of lo"'^ atoms 

-2 -1 
cm sec 
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Impact of lo plasma-torus electrons with the oxygen 
atoms determines both the lifetime of the oxygen atoms in 
the Jovian environment and the brightness of the oxygen 

O 

atom emissions in the 6300A lines. 

A two-dimensional distribution for the lo plasma-torus 
electrons was used in the calculation, and reflects the 
Voyager 1 encounter conditions as determined oy the UVS and 
plasma measurements. 

The two locations of Brown's observing aperature are 
indicated by the shaded rectangles , and his average measured 
value of 8 ± 4 Rayleighs is in agreement with the model 

calculation, if a slightly reduced satellite escape flux of 

9 -2-1 

5x10 oxygen atoms cm sec is adopted. 

27 -1 

Slide 2 This produces an ion creation rate of 2x10 ions sec 

having an initial spatial distribution relative to lo, shown 

in the next slide, as viewed from above the satellite plane. 

Note that the ion creation rate is concentrated near 

lo and mostly outside of its orbit, and also ahead of lo and 

m.ostly just inside its orbit. 

The outer contour corresponds to an ion creation value 
5 -2-1 

of 4.4 X 10 ions cm sec while the maximum value at lo 

7 -2 -1 

corresponds to an ion creation rate of 1.9x10 ions cm sec 
Assuming that equal numbers of sulfur ions would also 
be produced in a similar manner near lo, and that a plasma 
instability mechanism could thermalize the newly created co- 
rotational energies of these ions in about an hour, one would 
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Slide 3 


Slide 4 


Slide 5 


produce a hot electron source located just ahead of lo's 
orbital position with an energy input of edaout 2.6x10^^ 
watts or about 1 4% of the total energy radiated in the UV 
torus . 

This may provide an explanation for the lo-correlated 
energy source for the plasma torus recently discovered by 
Sandel from analysis of Voyager UVS data. 

In the next slide, the radial distribution of this 
ion-creation-rate is shown and indicates that ionization of 
oxygen atoms in the lo plasma torus is unsymmetric about lo's 
orbital location. 

This occurs because the oxygen atom lifetime is radial 
asymmetry about lo on the centrifugal equation of the plasma 
torus, as indicated in the next slide. 

The oscillation of the lo plasma torus about the satel- 
lite plane, which is not included in the present model cal- 
culations, will vary the oxygen atom lifetime at lo position 
by about a factor of 8. 

This oscillation will also vary the 6300A volume emission 
rate by a factor of about 6, so that the increased lifetime 
is approximately offset by a decrease in emission rate. 

Future model calculations, including the oscillating 
plasma torus, will refine the present model results. 

Titan Hydrogen Torus 

A Pre-Voyager model of the Titan hydrogen torus based 
upon analysis of Lyman-a emissions measured by the Earth- 
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orbiting Copernicus satellite and by the UV instrument of 
Pioneer 11 is shown on the next slide. 

A spatially uniform charge-exchange lifetime of 5600 hr 

7 9 

or 2x10 sec and a radial ejection flux of 1-3x10 H-atoms 

-2 - 1 

cm sec from Titan's exosphere with a mean emission speed 
of 2.0 km/sec were assumed to calculate the Lyman-a intensity 
pattern of a few hundred Rayleighs shown. 

The model is valid when there are no collisions between 
cloud atoms. 

The torus is nearly symmetric about Saturn except near 
Titan. 

The total height of the torus is about 6 Saturn radii 
and it has an apparent radial extent from about 8 to 28 
Saturn radii. 

UV observations of Voyager 1 measured a similar radial 
confinement of the Titan torus, but gave no indication of its 
vertical extent. 

The corresponding calculated radial distribution of the 

number of hydrogen atoms lost per 0.1 Saturn Radii per sec 

because of change exchange is given on the next slide and 

27 

represents a maximum ion exchange rate of about 2x10 ions 
-1 

sec 

The ion exchange rate is peaked at Titan's orbit and 
extends radially from 12 to 36 Saturn radii. A net ion 
creation rate of order 1% of these indicated value would be 
expected from electron impact ionization and photoionization. 

Recent Voyager 2 measurements of the Titan-hydrogen 
torus determined its total vertical extent to be about 
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12-14 Saturn Radii, about a factor of two larger than the 
pre-Voyager model results. 

Comparison of our pre-Voyager model parameters for the 
Titan torus with new values deduced from preliminary analysis 
of Voyager 2 UVS and plasma data are indicated on the next 
viewgraph. 

The major revision of the earlier model parameters 
results from a H-Atom lifetime 5 times larger. 

This, however, is offset by an increased torus volume 

of comparable size so that the effective ion-exchange rates 

27 

remain about the same, having a value of about 1x10 ions 
sec ^ . 

In addition to the improved values of the model para- 
meters, future models for the Titan hydrogen torus should also 
include the spatially non-uniform lifetime of oxygen atoms 
introduced inside of 8 Saturn Radii by the E-ring plasma 
torus and introduced outside of Titan's orbit by the solar- 
wind-magnetospheric boundary. 

This is illustrated in the last viewgraph, where the 
cross-sectional envelope of the Titan hydrogen torus is 
shown for atom escape speeds of 1.0, 1.5 and 2.0 km/sec. 

Furthermore, the validity of the model assumption of 
no collisions between cloud atoms may be called in question 
if the flux of H 2 or other molecules emitted by Titan is 

significantly larger than the deduced hydrogen atom flux of 

9 -2 - 1 

order 1x10 cm sec 
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ABSTRACT 


A new model has been developed for the sodium cloud of 
lo to describe its geocentric phase-angle dependent interac- 
tion with the solar radiation pressure arising from resonance 
scattering of sunlight in the and D, lines. Solar radia- 
tion pressure, acting along the Sun- Jupiter line, introduces 
an asymmetrically directed force into the otherwise cylindri- 
cal ly-symmetric gravitational three-body problem as described 
in a coordinate frame moving with the orbital motion of the 
satellite about Jupiter. The new model was developed to more 
quantitatively explore the hypothesis, originally presented 
in Paper I (Smyth, 1979) , that solar radiation pressure pro- 
vides a mechanism for explaining the east-west asymmetry in 
the cloud shape discovered by Goldberg et al . (1978) and the 
east-west cloud intensity asymmetry discovered by Bergstralh 
et al. (1975). From more recent observational data (Goldberg 
et al., 1980; Goldberg, 1981) these separately reported asym- 
metries now appear to be interrelated. Model calculations 
presented here confirm the original hypothesis and uncover 
explicit mechanisms for the observed east-west asymmetries. 
Model results presented for the cloud intensity do not, how- 
ever, include the spatially and time uependent sodium lifetime 
produced by the oscillation of the lo plasma torus about the 
orbital plane of the satellite. Modifications to these re- 
sults anticipated upon inclusion of the lo plasma torus in 
the model are, however, discussed and represent the only 
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major improvement to this model required before proper inver- 
sion of the east-west asymmetry data may be undertaken. 
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1. Introduction 


In an earlier paper (Smyth, 1979; hereinafter referred 
to as Paper I) , an explanation of the observed east-west 
asymmetries of the lo sodium cloud, based upon the effects 
of solar radiation pressure, was presented. These asymmet- 
ries were discovered by Bergstralh et al. (1975, 1977) and 
Goldberg et al. (1978) when Earth-based observations of the 
D-line emissions of the sodium cloud were compared for dia- 
metrically opposite satellite geocentric phase angles east 
and west of Jupiter. In Paper I, the force experienced by 
sodium atoms as they resonantly scatter sunlight in the D- 
lines and the ability of this force to alter the cloud atom 
orbits were explored. Orbit calculations presented showed 
that this solar radiation pressure produces a significant 
east-west asymmetry in the cloud shape. First, it acts in 
a time-dependent fashion, compressing the cloud near eastern 
elongation and expanding the cloud near western elongation, 
thereby providing the correct behavior to explain the east- 
west intensity asymmetry observed by Bergstralh et al. (1975, 
1977) . Second it causes the sodium cloud to tilt closer to 
Jupiter when lo is near western elongation that when lo is 
near eastern elongation, thereby providing an explanation 
for the east-west cloud shape asymmetry observed by Goldberg 
et al. (1978). An illustration from Paper I showing these 
effects is given in Figure 1. 
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The major objective of Paper I was to demonstrate simply, 
using orbit calculations restricted to the satellite plane, 
that solar radiation pressure provides a mechanism which may 
qualitatively explain the observed east-west asymmetry fea- 
tures. The object of this second paper is to further test 
this mechanism by providing quantitative and more realistic 
three-dimensional model calculations for the sodium cloud. 

To accomplish this, a sodium cloud model, including the in- 
herent time-dependent effects of solar radiation pressure 
on the orbits of cloud atoms, has been developed. This model 
represents a substantial modification of the earlier sodium 
cloud model of Smyth and McElroy (1977, 1978) and is presented 
in Section 2. 

In addition to solar radiation pressure, the actual in- 
tensity distribution of the observed sodium cloud vdll depend 
upon the initial velocity dispersion of sodium ejected from lo 
and also upon the spatial and temporal variations of the sodium 
lifetime in the Jovian environment. Knowledge of this veloc- 
ity dispersion should be obtained from a careful inversion of 
the sodium data using the above described cloud model into 
which the lifetime information and possible ion impact phen- 
omena (Brown and Schneider, 1982) have been properly incorpor- 
ated. The emphasis of the three-dimensional modeling presented 
in this paper will not, however, be data inversion- It is 
rather to provide the preliminary and necessary step of docu- 
menting the physical changes that occur in the sodium cloud 
intensity because of solar radiation pressure. 
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In this documentation, only the effects of solar radia- 
tion pressure on the lower velocity components (0-3 km sec”^) 
of the actual initial velocity dispersion of the ejected so- 
dium atoms are investigated. These components provide the 
dominant contribution to the intensity of the near lo sodium 
cloud. Higher emission velocities (3-15 km sec~^) that do 
not contribute dominantly to the intensity of the near lo 
sodinm cloud on the sky-plane (Smyth and McElroy, 1978) are 
not considered here. These higher emission velocities are, 
however, important, for example, in determining the asymmetric 
wings of the sodium line profile shape (Trafton, 1975; Trauger 
et al., 1976; Carlson et al., 1978; Macy and Trafton, 1980), 
in understanding the peculiar directional features seen by 
Pilcher (1980) and Goldberg et al. (1980), and in providing 
a source mechanism for atoms observed in the remote lo sodium 
cloud (Brown and Schneider, 1982) . In this documentation, 
changes in the cloud density, introduced because sodium atoms 
are lost in a spatially non-uniform manner by either ion 
impact or ionization through electron impact with the lo 
plasma torus, have also not been addressed. These changes 
represent a further numerical tailoring of the cloud density 
and D-line intensities within the three-dimensional sodium 
envelope which has already been altered by the action of 
solar radiation pressure. It is the altered intensity dis- 
tribution within this sodium envelope, without this additional 
numerical tailoring, that is the subject of this paper. The 
additional effects of this tailoring ar? ■ -rrently under 
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evaluation and will be considered in a future publication. 

A brief discussion of the anticipated effects of this tail- 
oring is, however, included in the following section. 

The discussion of model results for the east-west asym- 
metries is presented in Section 3 and Section 4. The discus- 
sion is divided into two sections for convenience of presen- 
tation, not because the results represent two unrelated phen- 
omena. In Section 3, model results for Goldberg's cloud 
shape asymmetry are discussed, while in Section 4, model 
results for Bergstralh's intensity asymmetry are discussed. 
Concluding remarks are presented in Section 5. 


2. Sodium Cloud Model 


An improved model for the lo sodium cloud is presented 
for calculation of its density and solar resonance scattered 
D-line intensities. The model is based upon following the 
trajectories of many orbits of atoms emitted by the satellite, 
subject to the combined gravitational fields of Jupiter and 
lo and the acceleration produced by resonance scattering in 
the D-lines. This involves numerically solving the equations 
of motion of each ^itom in thn e dimensions, for prescribed 
initial conditions, and including proper weight factors along 
each orbit to simulate ionization, elastic collisional less 
processes, and excitation phenomena. A description of the 
equations of motion and the resulting sodium cloud model, 
when the acceleration produced by solar resonance scattering 
is not included, was previously considered in detaiJ by Smyth 
and McElroy (1977, 1978). Similar approaches have also been 
adopted in other sodium cloud modeling efforts (Matson et al., 
1978; Goldberg et al., 1980). In these earlier papers, the 
motion of sodium atoms in the lo-Jupiter system required 
numerical solution of the circular restricted three-body 
problem in three dimensions. The presf’nt account provides 
an extension of these earlier papers by explicitly including 
the acceleration of sodium atoms produced by solar resonance 
scattering in the D-lines. 

The trajectories of sodium atoms in the improved model 
are governed by what will be called t.he modified circular 


8 



OF POOR QUALiTY 


restricted three-body equations of motion, given as follows; 
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b = b i 

sun 


(9) 

( 10 ) 


where b is the acceleration experienced by sodium atoms be- 
cause of resonance scattering in the D-lines, expressed in 
dimensionless form. The normal circular restricted three- 
cody equations of motion are recovered when b is equal to 
zero (see Smyth and McElroy, 1977). 
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The coordinates (x,y,z) measure the displacement of the 


sodium atoms with respect to the center of mass of the Ju- 
piter-lo system in units of the constant separation distance 
between the planet and satellite. Velocities (x,y,z) , accel- 
erations (x,y/z) and time t are defined in dimensionless 
fashion using the reciprocal of the angular frequency for 
motion of the planet and satellite around their center of 
mass as a characteristic time. Satellite and planet masses 
are denoted by mi and m 2 respectively. 

The coordinates (x,y,z) are defined in a non-inertial 
frame which rotates in such a manner to maintain the planet 
and satellite at fixed positions on the x-axis as shown in 
Figure 2 . These rotating coordinates are related to 
dimensionless inertial frane coordinates (S;,r!/S)» with right 


handed iinit vectors (i^/i^,i^), by the transformation 
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( 11 ) 


The inertial frame/ also with origin at the center of mass 
of the Jupiter-Io system/ is positioned with its C’^xis normal 
to the satellite orbit plane and its the plane 

defined by the Earth-Jupiter line and the ^-axis. For an 

Earth observer the direction of oositive is awav from the 
Earth toward Jupiter when the tilt angle of the satellite plane 

/V 

viewed from the Earth is zero. The unit vector i is directed 

sun 

away from the sun alona the Sun-Jupiter line. 
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The acceleration b is given by 

b = (bj + bj )/ 


G(mj + m^) 

. 7^ . 


( 12 ) 


where b and b are the contributions to the acceleration 
1 2 

of a sodium atom by resonance scattering in the and D 2 
lines respectively, where G is the gravitational constant 
and a is the constant seuaration distance between 
the planet and satellite. The acceleration b is made 
nondimensional by the characteristic acceleration given in 
brackets, which is effectively the gravitational acceleration 
experienced by a sodium atom in circular orbit of radius a 
about Jupiter. The magnitude of the acceleration b, as 
discussed in detail in Paper I, is as large as 0.01 to 0.02 
and depends upon the radial velocity v of the sodium atom 
with respect to the sun in the following manner: 


where 


m^,c 1 


Yi(v) 


(13) 


Tre ^ \ 





.-2 


(14) 


Here — ^ is the momentum change experienced by a sodium 
o 

atom of mas'" m^ upon absorbing a photon of energy hv^ 
in the line (i=l) or line (i=2) cuid is the corres- 
ponding rate at which photons are cdDSorbed. In the ex- 
pression for J^, R is the Sun-Jupiter distance in astronom- 

th. 

ical units, f^ is the oscillator strength of the i line. 
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ttF , is the solar continuum- level photon energy flux between 

the D and D lines at R = 1, and y. (v) is the fraction of 
12 1 

this solar continuum flux available to the sodium atom as 
it is Doppler shifted out of the bottom of the solar D-line 
Fraunhofer absorption feature by the instantaneous radial 
velocity v of the sodium atom relative to the Sun. 

The velocity v is the sum of the radial component of the 
velocity of the center of mass of the Jupiter-Io system 
u^ and the velocity of the sodiiim atom Ug^Qjjjin the inertial 
frame with respect to the Sun: 


where 


^'^cm ‘‘‘ ^aton ^ * ^un 


(15) 


u 


atom 




m2 ) 


( ? /H iQ) 


(16) 


and where (^,n,;) is given in terms of (x,y/Z) and (x,y,z) 
by differentiation of (11) with respect to t; 
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(17) 


The presence of the solar radiation acceleration in 

the equations of notion (1) ~ (3) nakes the differential 
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equations explicitly time-dependent in two ways. First there 
is the obvious explicit time-dependence in equations (1) - (3) 
introduced because of the changing vector direction of b in 
the (x,y,z) frame. Second, there is the complex time-depen- 
dent modulation of the magnitude of the acceleration vector 
b as described by (13) - (17) . The value of b varies more 
than an order of magnitude and is strongly dependent upon the 
satellite phase angle through the first term in (17) and is 
also dependent upon the instantaneous velocity of the sodium 
atom (x,y,z) relative to lo through the second term of (17). 
These two sources of explicit time dependence severely compli- 
cate the solutions of the equations of motion (1) - (3) . The 
natural circular symmetry of the gravitational problem that 
exists in the (x,y,z) frame when b = 0 is destroyed, thus caus- 
ing each atom ejected from lo with a given set of initial con- 
ditions to have a different trajectory for each different 
initial location of the satellite on its circular orbit.- 
This breakdown of the natural circular symmetry of the grav- 
itational problem caused by the introduction of the acceler- 
ation b gives rise to east-west differences in the cloud atom 
orbits discussed in Paper I and provides a mechanism to ex- 
plain the east-west asymmetries observed in the sodium cloud. 

Model calculations of the density and intensity of the 
sodium cloud are based upon determining the locus and instan- 
taneous velocity of cloud atoms in space resulting from con- 
tinuous ejection of sodium from lo. This locus of atom posi- 
tions for each of the individual initial conditions used to 
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describe this ejection process cannot be obtained by simply 
integrating the equations (1) - (3) once over the effective 
lifetime of the atom in the Jovian environment for one orbital 
location of lo as was accomplished in earlier modeling efforts 
where b is zero. In the present case, as the satellite is 
moving around its orbit continuously emitting sodium, the 
atoms ejected at different satellite orbital locations exper- 
ience a different force and therefore have different trajec- 
tories even though their initial conditions relative to the 
satellite are identical. The determination of the desired 
locus of atom positions in the (x,y,z) coordinates and their 
corresponding velocities therefore must be constructed, for 
each initial condition, from the solutions of the equations 
of motion (1) - (3) for a succession of different locations 
of the satellite along its orbit. The length of the angular 
section over which these successive solutions are required 
is determined by the effective lifetime of the sodium atoms. 
This angular section must also be properly pre-posit ioned along 
the satellite orbit to correctly define the cloud history at 
the time of observation. The number of satellite locations 
chosen within a given angular sector will be determined by 
the spatial resolution of the cloud density and intensity 
required in the model calculation. 

The equations of motion (1) - (3) may be taken to define 
a set of coupled first order differential equations which de- 
termine the time evolution of (x , y , z , x , y , z ) . These equations 
are solved numerically using a fourth order Runge-Kutta method. 
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The equations (1) - (3) exhibit no integrals of motion. 

Accuracy of the numerical solution is achieved by using 
carefully constructed se^<? of small tim ; steps for inte- 
grating each atom orbit. In performing these integrations, 
properly constructed weight factors along each orbit to 
simulate ionization, collisional loss processes, and D-line 
excitations of the sodium cloud atoms may also be calculated 
and used to compute the column density, D-line intensities, 
and D-line emission profiles of the cloud. 

In the model, the continuous ejection of sodium by lo 
is described by many orbits of atoms which are initially emitted 
from a spherical exobase. In the model calculations to fol- 
low, the atoms are initially emitted radially from the exobase 
taken to lie at an altitude of 780 km above the satellite 
surface. Other non-radial emission cases presently incorpor- 
ated in the model, such as ejection of sodium atoms from the 
satellite exobase by a Jupiter magnetospheric wind produced 
by the relative motion of the lo plasma torus past the satel- 
lite, are not considered here. These non-radial emission 
cases are of importance in studying data describing the asym- 
metric line profile shape of the sodium D-lines (Trafton, 

1975; Trauger, Roesler and Much, 1976; Trafton and Macy, 1977; 
Carlson et al., 1978), data describing the peculiar direc- 
tional features of the sodium cloud (Hartline, 1980; Pilcher, 
1980; Goldberg et al., 1980) and data describing the elastic 
ion-neutral collisions near lo that appear to populate the 
remote lo sodium cloud (Brown and Schneider, 1982). For 
model calculations considered in this paper, the exosphere 
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has 1298 source points distributed uniformly over the spher- 
ical surface. Radial emission of atoms from this exobase 
will, however, be restricted to monoenergetic ejection of 
sodium from the inner hemisphere of the satellite which 
involves only 685 of these source points. The inner hem- 
isphere is selected, since earlier modeling efforts (Smyth 
and McElroy, 1978; Matson et al., 1978; Macy and Trafton, 

1980) have shown it to be the dominant emission hemisphere 

to explain the two-dimensional intensity features of the sodium 
cloud. This assumption of inner hemisphere emission may 
actually reflect the presence of an extreme spatial asymmetry 
in the lifetime sink for sodium in the lo plasma torus to be 
discussed later. 

Ionizations of the sodium atoms along their orbits in 
the following model calculations are limited to the simplest 
description of a constant cutoff lifetime of 20 hours and no 
attempt is made, as discussed earlier, to address the more 
complex spatially non-uniform and time-dependent tailoring 
of the sodium cloud produced by its interaction with the lo 
plasma torus. The lifetime of sodium atoms in the lo plasma 
torus, based upon Voyager 1 data for the plasma density 
(Bridge, Sullivan and Bagenal, 1980; Bagenal and Sullivan, 

1981) and for the electron temperature deduced from UVS 
(Shemansky, 1980) and ^n situ measurements (Scudder, Sittler 
and Bridge, 1981) , is presented in Figure 3 and seen to be 
as small as about one hour near the satellite. This small 
lifetime value for sodium atoms is in marked contrast to the 
lifetime estimates of 15 to 20 hours (Smyth and McElroy, 1978) 
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and 28 hours (Matson et al., 1978) obtained from earlier mod- 
eling efforts and to the 20 hour value to be adopted here. 

In these earlier modeling studies, where no spatial varia- 
tion of the sodium lifetime was assumed, the value of the 
lifetime was essentially determined by properly modeling 
the spatial extent of the forward sodium cloud observed in 
Earth-based measurements. 

With the presence of the plasma torus since established, 
it is now clear that the forward extent of the sodium cloud 
is determined by ionization of these cloud atoms upon their 
secondary encounter with the plasma torus, somewhat ahead of 
lo, which occurs for flight times of order 20 hours (as illus- 
trated in Figure 4) . During the first encounter of the cloud 
atoms with the plasma torus (i.e., as they escape lo initially) 
the abundance of sodium atoms in the cloud is reduced in a 
manner which depends critically upon the time history of the 
atom orbits and their location in the plasma torus, which 
oscillates about the satellite plane. In Figure 3, the sod- 
ium lifetime at lo, for example, varies by a factor of three 
during this oscillation. This spatial and time-dependent 
destruction of the sodium density within the 20 hour cloud 
envelope will not be treated in this paper. Thi value of 20 
hours for the sodium lifetime is adopted, however, in order 
to preserve the correct dimensions for the forward sodium 
cloud. 

From Figure 4, it is relatively easy, however, to antic- 
ipate four effects of including the spatially non-uniform 
ionization of the lo plasma torus on the model calculations 
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presented here. First, note that the forward cloud, as 
determined by the lo plasma torus, will be significantly 
more elongated for western phase angles than for eastern 
phase angles, when cloud images are compared for diamet- 
rically opposite satellite locations about Jupiter. This 
more elongated behavior of the western cloud has been con- 
firmed in recent measurements of Goldberg et al. (1980) and 
Goldberg (1981) . The longer forward cloud in the west and 
the shorter forward cloud in the east will not, however, be 
a steady state phenomena, but will be seen as a time-depen- 
dent modulation about some mean state with an effective 
period of about 25 lo rotations or 44.25 days. This longer 
effective period results because of the non-resonance in the 
13 hour period of oscillation (with respect to lo) of the 
plasma torus about the satellite plane, and the 42.5 hour 
period of the asymmetric behavior of the sodium cloud intro- 
duced by the solar radietion pressure. 

Second, note that the density and D-line intensity of 
the diametrically opposite sodium clouds in Figure 4 will 
be modified differently by the lo plasma torus. The actual 
flight times for those sodium atoms which travel along the 
innermost edge of the forward east and west clouds, attaining 
a radial distance of 5.5 Jupiter radii at their second encoun- 
ters with the plasma torus, are about 28 and 24.8 hours 
respectively. This difference will produce a relative density 
enhancement in this portion of the forward east cloud. On the 
other hand, the outer edge (closest to the satellite orbit) 
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of the forward cloud in the east will experience a relative 
density reduction because of its more complete overlap with 
the lo plasma torus. The trailing portion of the east cloud 
outside of the satellite orbit will, however, have a relative 
density enhancement. This will occur because solar radiation 
pressure causes the atoms in this portion of the east cloud 
to move radially out of the plasma torus more rapidly than 
the corresponding portion of the west cloud, as can be seen 
in Figure 4. The net effect may well be an overall relative 
density enhancement for the east cloud, but further model 
computations are clearly required. These density changes 
will, of course, be time-dependent and may likely cause the 
mean east-west cloud shape asymmetry calculated in this paper 
to be modified in such a way as to affect the predicted value 
of the east critical phase angle more than the predicted 
value of the west critical phase angle. The critical phase 
angles are the satellite phase angles for which the sodium 
cloud on the sky plane appears to be distributed approximately 
symmetrically about lo. 

Third, note that because of the first ioniza,tion encoun- 
ter of the forward cloud atoms with the lo plasma torus, the 
ratio of the sodium density and D-line intensities very near 
lo to that in the forward portion of the cloud will be larger 
than in the present calculation. This will cause any east- 
west intensity asymmetry effects that originate near lo to be 
enhanced and will also require that a larger flux of sodium 
atoms be emitted from lo in order to properly populate the 
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elongated portion of the forward cloud. Brown (1981) has 
estimated that this increased sodium flux may be as large as 

c -ir26 ^ -1 

6 X 10 atoms sec 

Fourth, note that the extreme radial asymmetry of the 
sodium lifetime shown in Figure 3 and Figure 4 will provide 
a larger survival lifetime for sodium which moves inside the 
satellite orbital radius and populates the forward cloud, 
and a much smaller survival lifetime for sodium which moves 
outside the satellite orbital rad-Lus and populates the trail- 
ing cloud. This suggests that the apparent non-uniformity 
in the atom emission flux of sodium from lo (i.e., ejection 
from the inner hemisphere of lo with a spatially uniform 
lifetime of about 20 hours) deduced from earlier studies 
(Smyth and McElroy, 1978; Matson et al., 1978; Macy and 
Traftor. , 1980) might actually be explained by a symmetric 
ejection flux from lo with the very asymmetric spatial sink 
of the plasma torus. Comparison of two-dimensional sodium 
cloud images with very preliminary results calculated with 
our model including a non-oscillating plasma torus suggests, 
however, that the asymmetric sink is not strong enough. In- 
clusion of the oscillating motion of the plasma torus in the 
model may, however, provide the necessary time-average en- 
hancer..ent of the sink asymmetry needed to allow calculations 
assuming symmetric emission from lo to be consistent with the 
observed sodium images. 

Initial velocities of atoms emitted from the satellite 
exobase in the following model calculations are limited to 
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values less than 3 km sec ^ for two reasons. First, the 
mean emission velocity associated with the formation of the 
forward portion of the near lo cloud, for the brightness 
levels greater than about 1000 Rayleighs viewed by Goldberg 
et al. (1978), is 2.6 km secT^ (Smyth and McElroy, 1978). 

This value of the mean emission velocity determines the 
inclination angle of the forward cloud, the angle between 
the central axis of the forward cloud and the tangent line to 
the satellite orbit at lo (see Figure 1) which was shown in 
Paper I to be modulated by solar radiation pressure. Larger 
mean emission velocities produce mean inclination angles that 
are too large whereas smaller mean emission velocities pro- 
duce mean inclination angles that are too small. Second, 
sodium atoms ejected from the exobase, with emission veloci- 
ties slightly in excess of the satellite escape speed at that 
altitude, provide the dominant contribution to the D-line 
intensity in the spatial region close to lo observed by 
Bergstralh et al. (1975, 1977). This behavior is illustrated 
in Figure 5 where the relative D 2 intensity contribution to 
their 3x8 arc sec slit centered on the satellite is shown 
for monoenergetic model calculations, excluding solar radia- 
tion pressure, with emission velocities ranging from 1.8 km 
sec ^ to 10 km sec The low velocity of 2.0 km sec ^ 
dominates since atoms with this emission velocity move very 
slowly with respect to lo upon escape. Lower emission 
velocities have a larger ballistic component and contribute 
less, whereas higher emission velocities move much more 
rapidly away from lo and contribute less. 
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3. Model Results: Goldberg's Cloud Shepe Asymmetry 

The sodium cloud has been observed by Goldberg et al. 
(1978) to exhibit an east-west asymmetry in the spatial 
distribution of the D-^ine intensities (i.e., in the cloud 
shape) with respect to the geocentric orbital phase angle 
of the satellite. Their observations showed that the ob- 
served intensity pattern of the cloud, which is its projec- 
tion on the sky plane, did not produce a mirror image for 
satellite orbital phase angles separated by 180 degrees. 

To initially quantify and characterize one of the differ- 
ences observed in comparison of these east and west cloud 
images, the critical phase angles were selected as param- 
eters. The cloud exhibits two critical phase angles which 
occur for those satellite phase angles where the sodium 
cloud on the s)cy plane appears nearly symmetric about lo as 
the satellite approaches either the eastern or western elong- 
ation points of its orbital (i.e., 90 degrees or 270 degrees, 
respectively) . 

The first estimates of the eastern and western critical 
phase angles were originally reported by Goldberg et al. 

!1978) as 65 degrees and 230 degrees respectively. These 
values were defined quantitatively by measuring, from their 
two-dimensional cloud images, the integrated intensities seen 
in two simulated slits located 9 arc sec on either side of 
lo, and subsequently by determining that satellite phase 
angle for which the ratio of the two slit intensities was 


22 



unity. This method led to a rather well defined value for 
the critical phase angle in the west, but to a rather ill 
defined value in the east due to scatter. The less distinct 
behavior of the data in the east seemed to reflect a funda- 
mental east-west difference in the sodium cloud, not a dif- 
ference in the quality of the east and west da^a. Later 
estimates made with additional observations (Goldberg, 1979) 
indicated chat the critical eastern phase angle might even 
be as large as 70 to 75 degrees. This then suggested that 
in the east the central axis of the forward sodium cloud 
was inclined only abou’j 15 to 20 degrees with respect to the 
line drawn tangent to the satellite orbit at lo, whereas in 
the west this inclination angle was about 40 degrees. The 
sodium cloud would then nave an east-west phase lag asymmetry 
of about 20 to 25 degrees. 

The accuracy of these estimates for the east and west 
critical phase angles and hence for the east-west phase lag 
asynuTietry has been, however, somewhat difficult to determine 
with exactness from the data of Goldberg et al. (1978) and 
Goldberg (1979) , since changes in the apparent cloud geometry 
during their 2 to 3 hour long integrations correspond to 
orbital angles of 17 to 25 degrees and reduce the spatial 
resolntioo of their cloud images. To rectify this situation, 
Goldberg (1981) has recently made several new measurements 
using an improved instrument with an integration time of 
about 10 minutes. His preliminary analysis of these new 
measurements, based upon visual comparisons, indicates that 
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the western critical phase angle is distinct and occurs near 
235 degrees, whereas the eastern critical phase angle is 
significantly less distinct and occurs somewhere in the 
angular interval from 55 to 70 degrees. The exact value 
of the east-west phase lag asymmetry at present, therefore, 
remains undetermined and may or may not be as large as in- 
itially suggested. The one significant east-west asymmetric 
cloud feature that does occur both in these new measurements 
and in their older observations is the obvious presence of a 
distinct critical western phase angle and the noticeable 
visual absence of a distinct critical eastern phase angle. 
Quantitative analysis of existing data and acquisition of new 
data will be required to determine more accurately the critical 
eastern phase angle and to understand possible variations 
of its value introduced by both the oscillating lo plasma 
torus and non secular variations in the local plasma. Fur- 
ther analysis of these data should also be useful in iden- 
tifying and characterizing other salient features of the 
east-west cloud distibutiou asymmetry, such as the greater 
elongation of the forward portion of the sodium cloud ob- 
served for western phase angles by Goldberg et al. (1980) 
and Goldberg (1981) . To provide a step forward, based upon 
our current understanding, modeling analysis in this section 
will be restricted to exploring and documenting the east- 
west phase lag asymmetry introduced by the solar radiation 
pressure, excluding the spatially non-uniform and time- 
dependent ionization of sodium atoms by the lo plasma torus. 
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The preliminary results of calculations in the satellite 
plane, originally presented in Paper I and repeated here in 
Figure 1, demonstrate that solar radiation pressure in the 
D- lines causes the envelope of the sodium cloud to behave in 
a manner similar to the observations of Goldberg et al. (1978) . 
Changes in the cloud shape in Figure 1 suggest an east-west 
phase lag of about 15 degrees. It is clear, however, that 
such an estimate is only preliminary and that the actual 
three-dimensional distribution of sodiiim atoms within the 
cloud envelope must be calculated in order to determine a 
more accurate phase lag angle. Suitable three-dimensional 
model calculations for this purpose are presented below. 

Three-dimensional model calculations, illustrating the 
effects of solar radiation acceleration on the D 2 intensity 
distribution of the sodium cloud viewed parallel to the sat- 
ellite plane, are given tor satellite phase angles of 55 and 
235 degrees in Figure 6 and Figure 7, respectively. In these 
model calculations, which may be compared directly with the 
simpler results of Figure 1, sodium was emitted radially from 
the inner hemisphere of lo's exobase (2600 km radius) with a 
velocity of 2.6 km sec ^ and with a cutoff lifetime of 20 

hours in the Jovian environment. Each outer contour value 

8 “2 “1 

is 916 Rayleighs for a surface flux of 10 atoms cm sec , 
with contours spaced by 367 Rayleighs, increasing inward 
toward lo. The results of Figure 6 and Figure 7 show that 
changes in the D 2 intensity pattern of the sodium cloud are 
not very sensitive to the illumination angle of the sun rel- 
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ative to the Earth- Jupiter line of sight. The sun angle 
(i.e., its phase angle plus 180 degrees) is con- 
fined to values between about 168 and 192 degrees as the 
earth moves about the Sun. For a sun angle of 180 degrees, 
the Earth-Jupiter line of sight and the Sun-Jupiter axis of 
Figure 1 are identical. For a satellite phase angle of 55 
degrees, note that in Figure 6 the effect of solar radiation 
is to rotate the central axis of the forward elongated sodium 
cloud further away from Jupiter as anticipated in the results 
of Figure 1. With solar radiation acceleration included in 
the model calculation, the inner edge of the cloud (right 
portion) in Figure 6 is therefore closer to lo while the 
outer edge (left portion) is extended further from lo. For 
a phase angle of 235 degrees, note that in Figure 7 the solar 
radiation pressure has effectively rotated the sodium cloud 
axis more toward Jupiter, as anticipated in Figure 1. With 
solar radiation acceleration included, the outer edge of the 
cloud (right portion) in Figure 7 is therefore closer to lo 
while the inner edge of the cloud (left portion) clearly re- 
veals that the axis of the forward elongated sodium cloud 
has already begun to swing through our line of sight. 

The model calculations of Figure 6 and Figure 7 indicate 
that solar radiation pressure introduces an east-west phase 
lag asymmetry into the spatial distribution of the sodium 
cloud. This can be illustrated K.ore clearly as follows. In 
the absence of solar radiation acceleration, model calcula- 
tions of the D 2 intensity pattern of the sodium cloud for 
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diametrically opposite satellite phase angles are presented 
in Figure 8, where the sun angle is 180 degrees and the 
other model parameters and contour levels are the same as 
in Figure 6 and Figure 7. Note that the east and west in- 
tensity patterns are effectively mirror images of each other, 
and that tne axis of the forward elongated cloud can be seen 
to swing through our line of sight at a phase angle of about 
55 degrees in the east and 235 degrees in the west, yielding 
no phase lag. The observed changes in the intensity pattern 
of the cloud with satellite phase angle are almost completely 
due to the changing geometric viewing perspective of the 
cloud projected onto the plane of the sky. In Figure 9, the 
D 2 intensity patterns of the sodixam cloud, where the effects 
of solar radiation acceleration have been included, are pre- 
sented for the identical diametrically opposite satellite 
phase angles, model parameters and contour levels used in 
Figure 8. Near eastern elongation, note that the outer con- 
tour of the sodium cloud appears symmetric about lo near 
a satellite phase angle of between 60 and 65 degrees, where- 
as near western elongation the cloud appears symmetric at a 
phase angle of 230 degrees. This produces an east-west phase 
lag asymmetry between 10 and 15 degrees. If, on the other 
hand, the distribution of the D 2 intensity within the outer 
contour in Figure 9 is used as a criterion, approximate 
symmetry occurs for a chase angle of 65 degrees in the east 
and 235 degrees in the west. This produces an east-west 
phase lag asymmetry that is also about 10 degrees. Note 
tha;. in Figure 9, the passage of the forward cloud axis 
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through our field of view is much le^^s less distinct in 
the east than in the west. This makes it more difficult 
to accurately determine the eastern critical satellite phase 
angle, a persistent characteristic noted in the observational 
data of Goldberg et al. (1978, 1980) and Goldberg (1979, 1981). 

The results in Figure 8 and Figure 9 are calculated under 
the assumption of a monoenergetic emissi'^n velocity of 2.6 
km sec Results for lower emission velocity components 

that would be included if a more realistic initial emission 
velocity distribution were adopted have also been calculated. 
For emission velocities of 2.0 km sec 2.2 km sec ^ and 2.4 
km sec the intensity patterns on the sky plane are again 
essentially mirror images, in the absence of solar radiation 
pressure, when the east and west clouds are compared for di- 
ametrically opposite satellite phase angles. The critical 
phase angles increase with decreasing emission velocity, 
having a value, in the absence of solar radiation pressure, 
in the east between 55 and 60 degrees for an initial veloc- 
ity of 2.4 km sec ^ and a value of about 60 degrees for in- 
itial velocities of both 2.2 km sec ^ and 2.0 km sec 
With solar radiation pressure included, these model calcu- 
lations have the same general east-west character as the 
results presented in Figure 9. Based upon an outer contour 
criterion, the critical eastern phase angle is about 65 de- 
grees for the emission velocities of 2.4 km sec ^ and below, 
while the critical western phase angle is between 230 and 
235 degrees. This provides an east-west phase lag asymmetry 



28 


of about 10-15 degrees. Based upon a criterion for r 
symmetric distribution of intensity within the outer contour, 
the critical eastern phase angle is about 70 degre*^-? while 
the critical western phase angle is about 235 degrees for 
emission velocities of 2.4 km sec ^ and lower, producing a 
15 degree east-west phase lag asymmetry. 

Model results for emission velocities of 2.2 km sec 
2.4 km sec ^ and 2.6 km sec ^ are compared in Figure 10 and 
Figure 11 for diametrically opposite phase angles east and 
west of Jupiter These D 2 intensity images include the 
effects of solar radiation pressure and are calculated for 
the same contour levels, emission conditions, and lifetime 
conditions chosen in Figures 6-9. As the emission velocity 
decreases, the cloud is more closely confined near lo and 
provides a rapidly increasing intensity near the satellite. 
The exact east and west critical phase angles and correspond- 
ing east-west phase lag asymmetry predicted from these model 
calculations will therefore depend to some extent upon the 
shape of the initial velocity distribution. From the above 
calculations, however, it is clear that a phase lag of about 
10 to 15 degrees would be present even if an initial velocity 
distribution were adopted. For such an adopted distribution, 
the passage of the cloud through its eastern critical phase 
angle would also appear to occur slower than the passage of 
the cloud through its western critical phase angle. The 
critical phase angle in the east should be within the angular 
range from 55 to 70 degrees, while in the west it should be 
between 230 and 235. 
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The presence of a distinct critical western phase angle 
and of a significantly less distinct critical eastern phase 
angle in the observations of Goldberg (1981) suggests that 
the non-uniform ionization of sodium introduced by the lo 
plasma torus oscillating about the satellite plane may tailor 
and modify, but will not eliminate, the basic characteristics 
of the east-west cloud distribution asymmetry illustrated by 
the model calculations in this section. Further quantitative 
analysis of the data acquired in the 1980-1981 apparition by 
Goldberg (1981) is warranted. Such analysis should also pro- 
vide important information for modeling the alternating 
north-south asymmetry discovered by Trafton and Macy (1975) 
and discussed further by Trafton (1980) . 
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4. Model Results; Bergstralh's Intensity Asymmetry 

The sodium cloud has been observed by Bergstralh et al. 
(]975, 1977) to exhibit an east-west asymmetry in its abso- 
lute intensity. In their observations, the D-line intensi- 
ties of the sodium cloud were measured as a function of phase 
angle through a 3" x 8" slit centered on the satellite with 
the long dimension of the slit oriented approximately perpen- 
dicular to the projection of lo's orbital plane. This obser- 
vational geometry is ili-^ustrated in Figure 12. The long 
dimension of the slit is only slightly larger than the diam- 
eter of the Lagrange sphere of lo, within which the gravita- 
tional field of the satellite dominates the planetary gravi- 
tational field. The spatial region of the cloud sampled by 
the slit is therefore composed of a mixture of the dense 
sodium that likely exists within the Lagrange sphere and the 
somewhat less dense sodium that occurs, along the line of 
sight, in the foreground and background of this sphere. In 
this spatial region, sodium atoms emitted from the lo exobase 
with velocities very near the satellite escape speed are ex- 
pected to provide the dominant contribution to the slit in- 
tensity as shown earlier in Section 2 and in Figure 5. 

The D-line intensities of the cloud, measured through 
this observing slit, were shown by Bergstralh and colleagues 
to be correlated with lo's orbital position around Jupiter. 
This correlation made it possible for them to identify reson- 
ance scattering of sunlight by sodium atoms as the dominant 
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excitation mechanism for the cloud. In addition, the measured 
intensity of the cloud with the satellite east of Jupiter was 
shown, in their initial paper (Bergstralh, Matson, and John- 
son, 1975) to be consistently about 50% greater than the mea- 
sured intensity with the satellite west of Jupiter. In their 
second paper (Bergstralh et al., 1977), which included addi- 
tional observations and a presentation of more select data, 
this east-west intensity asymmetry was more accurately estim- 
ated to be 20 to 25%. 

In Paper I, this intensity asymmetry was associated with 
the ability of solar radiation pressure to compress the east 
cloud and expand the west cloud envelope, as illustrated in 
Figure 1. To investigate more fully how solar radiation pres- 
sure provides a mechanism to explain this intensity asymmetry, 
a number of three-dimensional model calculations have been per- 
formed and are presented here. The and D 2 intensities seen 
through a 3" x 8" slit centered on lo have been calculated as 
a function of both the satellite phase angle and the initial 
speed of sodium atoms emitted from the inner hemisphere of 
lo' s exobase. In all of these model computations, the two- 
dimensional distributions of both the column density and the 
intensity of the cloud on the sky plai as well as the shape 
of the sodium D-line profiles, were simultaneously calculated, 
both excluding and including solar radiation pressure, to 
understand more fully the basic phenomena. 
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Model calculations for the relative D 2 intensity seen 

through the 3" x 8” viewing slit, taken here to have nominal 

3 4 

dimensions of 9 x 10 km by 2,4x10 3cm, are presented in 
Figure 13 and compared with the observations of Bergstralh 
et al. (1975). The relative D 2 intensity of the model cal- 
culations is shown as a function of the satellite phase angle 
for four different monoenergetic emission velocities spanning 
the range from 2.0 1cm sec ^ to 2.6 km sec In each mono- 

energetic model calculation , the same number of atoms were 
emitted from lo's exobase. These relative intensities were 
calculated for sodium atoms emitted from the inner hemisphere 
of the exobase and include the shadow effect of the disk of 
lo blocking the cloua behind it. The tilt angle of the sat- 
ellite plane is assumed to be zero relative to the observer. 
The sun, observer and Jupiter are assumed to be along the 
same straight line. Intensities calculated for a satellite- 
plane tilt angle of a few degrees differ only slightly. 

The model calculations of Figure 13 clearly show that 
solar radiation pres^ •’■^e introduces a large east-west inten- 
sity asymmetry for the initial emission velocity of 2.0 1cm 
sec ^ and a somewhat less pronounced asymmetry for the 2.1 
km sec ^ emission velocity. Only a small asymmetry 
is present for the higher emission velocities of 2.4 
km sec ^ and 2.6 km sec The maximum intensity of the 2.0 

1cm sec ^ results is 4,38 times larger than the 2.6 km sec ^ 
results in the east and 3.48 times larger in the west. The 
relative intensities for the different velocities in 


33 


Figure 13 have a maximum value near 60 degrees in the east 
and near 230 degrees in the west. In both the east and west, 
these peak intensities correspond to the enhancement of the 
sodium cloud density produced by the forward portion of the 
cloud as it passes through the line of sight of the viewing 
slit. The peak intensity is larger and more pronounced for 
the lower emission velocities since the forward cloud is then 
much more narrow and is much more tightly confined near lo. 
This spatial behavior of the forward cloud, in the absence 
of solar radiation pressure, was discussed earlier by Smyth 
and McElroy (1977) , who at that time noted a similarity be- 
tween their calculated column density peaks and the peaks in 
the data of Bergstralh et al. (1975) . This similarity is 
also present in the model results of Figure 13, although 
before an accurate comparison can be made with the observa- 
tions, the additional modifications introduced by the spa- 
tially non-uniform lifetime of the lo plasma torus must be 
included in the model calculations. 

The east-west ratio of the relative D 2 intensity for 
the 2.0 km sec ^ and 2.1 km sec ^ emission velocities is 
presented in Figure 14, where the results of including and 
excluding the effects of Io*s shadow are both indicated. 

For an emission velocity of 2.0 km sec ^ the ratio provides 
a minimum east-west asymmetry of 17 percent over most of the 
angular range when the shadow is included, and at a 70 degree 
phase ang]e the ratio has a peak value of over 1.8. The 2.1 
km sec ^ ratio is less pronounced, but still provides a sig- 


niJicant east-west asymmetry over much of the phase angle 
interval. The general qualitative behavior of the in- 
tensity ratio exhibited in Figure 14 can be simply under- 
stood in terms of the action of solar radiation pressure 
upon the orbits of sodium atoms. This qualitative explan- 
ation will now be discussed in some detail to clarify both 
the nature of the physical effects involved and the contri- 
bution of the phase-lag asymmetry phenomena to the intensity 
asymmetry. 

In Figure 14, the D 2 intensity ratio is greater than 
unity for east/west phase angles in the range from 0/180 to 
40/220 because the column density of the east cloud in the 
viewing slit is larger than that of the west cloud by approx- 
imately 20% for an emission velocity of 2.0 km sec The 

radial velocity v of most sodium cloud atoms with respect to 
the sun is larger, however, in the west cloud than in the 
east cloud, for this phase angle interval, so that the 
intensity ratio is reduced somewhat from this column density 
ratio through the factor y^(v) in equation (14). This larger 
velocity of sodium atoms in the west cloud is also evident in 
sodium line profiles when model results are compared for di- 
ametrically opposite satellite phase angles. For an east/ 
west phase angle near 50/230 in Figure 14, the D 2 intensity 
of the west cloud in the viewing aperature has reached its 
maximum value, because of the phase lag asymmetry alignment 
of the forward cloud. The east cloud, however, reaches its 
maximum D 2 intensity value of about 10 degrees later and 
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maintains intensity values near this maximum level over a 
longer phase angle interval, beyond its peak, than does the 
west cloud (see the model calculation in Figure 13) . This 
gives rise in Figure 14 to the dramatic increase in the D 2 
intensity ratio within the east/west phase angle range from 
50/230 to 120/300, and to the asymmetric slope of the D 2 
intensity ratio about its peak value located at a phase 
angle of 70/250. 

For larger east/west phase angles between 120/300 and 
165/345, the D 2 intensity ratio decreases and eventually 
becomes less than unity. In this angular range, most atoms 
now seen in the viewing slit are near lo and not in the for- 
ward cloud. Solar radiation pressure, in this angular range 
so acts on these cloud atoms to elevate the column density 
of the west cloud by as much as 18% for the 2.0 km sec ^ 
emission velocity and as much as 4“' for the 2.1 )cm sec ^ 
emission velocity. In the west cloud, the radial velocity 
V of most atoms with respect to the sun is also enhanced 
over the east cloud and, through the factor r^(v) in equa- 
tion (14), causes the D 2 intensity ratio to be reduced even 
further . 

The significant enhancement of the calculated column 
density of the east cloud relative to the west cloud within 
the phase angle range 0/180 to 40/220 in Figure 14, may be 
attributed to two factors. The first one is that solar radi 
ation pressure produces an east cloud that is spatially more 
tightly confined near lo and hence has elevated densities. 
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The second and more important factor is that solar radiation 
pressure modifies the escape process of atoms from lo so that 
the total number of atoms within the east cloud is larger. 
This modification of the escape process is also largely re- 
sponsible for the elevation of the column densities of the 
west cloud relative to the east cloud within the east/west 
phase angle range from 135/315 to 170/350. 

With regard to the first factor, the density enhance- 
ments of the east cloud and the density reductions of the 
west cloud within the phase angle range 0/180 to 40/220 in 
Figure 14 are more pronounced for atom emission velocities 
near 2.0 km sec This occurs because many of these atoms, 
shortly after emission from the exobase, experience a near 
balance between the effective gravitational field of Jupiter 
and lo. This near balance causes the relative magnitude of 
the solar radiation acceleration experienced by the sodium 
atom to be favorably amplified. The vector directio’^ of the 
solar radiation acceleration, being oppositely oriented in 
the east and west cloud, then asymmetrically and somewhat 
abruptly alters the orbits of the cloud atoms. Future evo- 
lution of these altered atom orbits then produces the density 
enhancement of th east cloud ana density reduction of the 
west cloud. Fot uhe higher emission velocities of 2.4 km 
sec ^ and 2.6 km sec this amplification effect is much 
smaller since che sodium atoms, moving more rapidly away 
tiom lo, sp3=nd less time near the Lagrange sphere. The 
resulting east-west density changes for these higher emis- 
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Sion velocities are therefore much less pronounced and result 
primarily from the smooth and continuous action of solar ra- 
diation pressure over the 20 hour lifetime of the atoms. In 
the case of the higher emission velocities, the compression 
cf the east cloud and the expansion of the west cloud, as 
shown in Figure 1, are, in fact, to first order compensated 
by the column integration along the line of sight and con- 
tribute little to the intensity asymmetry. 

With regard to the second and more important factor, 
the ability of solar radiation pressure to modulate the 
total number of sodium atoms within the cloud occurs essen- 
tially for those emission velocities near 2 km sec ^ where 
the escape of atoms is energetically possible only from a 
limited surface area of the satellite exobase. This ener- 
getically possible area of the exobase defines the window 
in the Lagrange sphere surrounding lo through which the 
atoms escape. It is this effective exospheric area, from 
which the emitted atoms are able to escape from the satel- 
lite, that is modulated by the action of solar radiation. 

The Lagrange window in lo's orbital plane for atoms 
emitted from the exobase of the satellite has been discussed, 
in the absence of solar radiation pressure, by Smyth and 
McElroy (1977). One of their iilustr ions is shown in Fig- 
ure 15 for radial emission velocities of 1.8 km sec~^, 2.0 km 
sec ^ and 3.0 km sec a 1.8 km sec ^ velocity, the 

Lagrange window is closed except near the canter of the 
inner satellite hemisphere (i.e., at 180°) and near the 
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center of the outer satellite hemisphere {i.e., at 0°). 

The window is preferentially open in this direction because 
it is the easiest direction for escape of atoms, since the 
two col linear Lagrange points, and L 2 , nearest to lo lie 
along the line connecting lo and Jupiter (see Figure 2) . 

For an emission velocity of 2.0 km sec the Lagrange win- 
dow is closed on a small sector of the leading inner quadrant 
and a small sector of the trailing outer quadrant of the exo- 
base. For a velocity of 3.0 km sec the Lagrange window 
is completely open, since all atoms emitted radially from 
the exobase now have sufficient energy to escape lo regard- 
less of their initial surface location. In the absence of 
solar radiation pressure, the percentage of atoms emitted 
from the inner hemisphere of the exobase that do not escape 
lo is determined from model calculations to be 60.3, 28.6, 9.6, 
1.8 and 0.19 for monoenergetic emission velocities of 1.8 km 
sec 1.9 km sec 2.0 km sec 2.1 km sec ^ and 2.2 km 
sec These percentage values are, of course, independent 

of the phase angles at which tiie satellite is observed. 

The ability of solar radiation pressure to modulate, 
as a function of the satellite phase angle, the percentage 
of atoms not escaping from lo, having been initially emitted 
radially from the inner hemisphere of the exobase, is illus- 
trated in Figure 16. Model results, wher j solar radiation 
pressure has been noth included and excluded, are compared 
for monoenergetic emission velocities of both 2.0 ’em sec ^ 
and 2.1 km sec The functional behavior exhibited is 
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essentially similar in nature for both emission velocities. 
Solar radiation pressure enhances the escape of sodium atoms 
for smaller satellite phase angles in the east (i.e., be- 
tween 0 and 68 degrees for 2.0 km sec ^ and between 0 and 
98 degrees for 2.1 km sec and for large phase angles in 
the west (i.e., between 208 and 360 degrees for 2.0 km sec ^ 
and between 291 and 360 degrees for 2.1 km sec ^) . This 
behavior accounts, in Figure 14, for the earlier noted en- 
hancement of the column density of the east cloud relative 
to the west cloud in the east/west phase angl>^ range of 
0/180 to 40/220, and of the west cloud relative to the east 
cloud in the east/west phase angle range of 135/315 to 170/ 
350. The general behavior exhibited by the effect of solar 
radiation pressure in Figure 16 may also be anticipated on 
geometrical grounds by simply considering, as a function of 
the observed phase angle of the satellite, the 20 hour time 
history of the orientation of the solar radiation accelera- 
tion vector and that portion of the inner hemisphere exobase 
froai which atom escape is modulated. 

The maximum observed D 2 intensity reported by Bergstralh 
et al. (1977) using their best data base was 241 k Rayleighs. 
This was measured for a satellite phase angle near eastern 
elongation and was determined by assuming that only a 3 x 3 arc 
second portion of their observing slit was uniformly illumi- 
nated. If the complete 3x8 arc second slit is considered, 
this uniform intensity is reduced to about 90 k Rayleighs. If 
this intensity were produced by ths maximum D 2 intensity given 
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by the 2.0 km sec ^ calculation of Figure 13, a satellite flux 
8 ~2 —1 

of order 1 x 10 atoms cm sec would be required, in agree- 
ment with an earlier estimate by Smyth and McElroy (1978) . 

8 "2 ~ 1 

This flux value of 1 x 10 cm sec will, however, be re- 
vised upward, as noted earlier by Brown (1981) , upon inclu- 
sion in the model calculation of the highly spatially, non- 
uniform ionization of sodium atoms by the lo plasma torus. 

Modification of the results of Figure 13 and Figure 14 
by including the presence of the lo plasma torus in the new 
model calculations may be anticipated on simply geometrical 
grounds. The forward sodium cloud for an emission velocity 
of 2.0 km sec ^ will be largely confined between the inner 
and outer radial boundary of the plasma torus, in contrast 
to the 2.6 km sec ^ results illustrated in Figure 4. The 
continual exposure of this sodium to the ionization of the 
plasma torus will cause the atom population in the forward 
cloud to be significantly diminished relative to the atom 
population near lo. The dramatic increase in the model 
calculated D 2 intensity ratio in Figure 14 for an east/west 
phase angle of 70/250, which occurs because of the passage 
of the forward cloud through t.he field of view of the ob- 
serving slit, will therefore be reduced by this ionization 
of the plasma torus. The extent of this reduction will 
depend upon the mean plasma conditions, whicr. have oeen 
observed to change over a time period of several months 
(Sandel et al., 1979), and also upon the time modulation 
of these mean conditions by the oscillation of the lo plasma 
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torus about the satellite plane. Inversion of the east-west 
asymmetry data of Bergstralh et al. (1975, 1977) using an 
improved model, which properly includes these plasma torus 
effects, may then be expected to lead to meaningful results. 
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5. Concluding Remarl. 

A three-dimensional model for lo's sodium cloud incor- 
porating the full effects of solar radiation acceleration 
associated with solar resonance scattering of atoms in the 
and D 2 lines has been developed. The model has been 
applied to interpret two different but related east-west 
asymmetries exhibited in the near lo sodium cloud in data 
obtained from ground-based telescopes. The major objective 
cf the modeling analysis presented was to document the 
physical changes produced in the sodium cloud by the solar 
radiation acceleration, as a function of the satellite phase 
angle, and to further clarify beyond the efforts of Paper I, 
how these changes provide an explanation of the observed 
east-west asymmetries. Model results presented did not 
include the effects of spaticilly non-uniform and time var- 
iable interactions of cloud atoms with the lo plasma torus. 
Anticipated modifications of these results by the playma 
torus were, however, discussed and are expected to introduce 
quantitative rather than qualitative changes. The develop- 
ment of an improved model to include the plasma torus in 
a quantitative fashion is currently in progress. 

Model results for the observed east-west cloud shape 
asymmetry of Goldberg et al. (1978, 1980) and Goldberg (1979, 
1981) wer€i presented in Section 3. These results, calculated 
for different monoenergetic emission velocities, reproduced 
the correct phase angle behavior of the observations, indi- 
cating a distinct critical phase angle in the west between 
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230 and 235 degrees and a significantly less distinct crit- 
ical phase angle in the east between 55 and 70 degrees. 

These results suggest an east-west phase lag asymmetry of 
10 to 15 degrees. The general behavior of the cloud shape 
asymmetry and the value of the phase lag asymmetry predicted 
by the model calculations are not critically dependent upon 
the exact distribution of the emission velocities between 

2 km sec ^ and 3 km sec ^ (the velocity range which contrib- 
utes most dominantly to the intensity of the near lo cloud) , 
because these effects depend primarily upon the simple geo- 
metrical orientation of lo's orbit and the solar radiation 
force vector. 

Model results for the observed east-west intensity 
asymmetry of Bergstralh et al. (1975, 1977) were presented 
in Section 4. Intensity contributions seen through their 
3x8 arc sec viewing slit centered on lo were shown to be 
do'-. ’Dated by atoms with emission velocities less than about 

3 Km sec with the most important contribution occurring 
for an emission velocity of 2 km sec a value very near 
the satellite escape speed at the exobase. Model results 
near 2.C km sec ^ pro\'ide a distinct east-west intensity 
asymm.etry over most of the satellite phase angle range and 
compare favorably v’ith the observational data. Solar radia- 
tion pressure produces this intensity asymmetry by modulat- 
ing, as a function of satellite phase angle, the size of the 
Lagrange window through which atoms are able to escape the 
gravitational field of lo. For atoms with emission veloci- 
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ties larger than about 2.2 km sec the perturbing effect 
of solar radiation pressure is not strong enough to modulate 
the size of the Lagrange window. East-west intensity asym- 
metry effects for these larger velocities are introduced 
primarily by the different overlap of the cloud envelope 
and the lo plasma torus. This overlap is time-dependent, 
having an overall effective period of 25 lo rotations or 
44.25 days, with the mean results likely favoring a rela- 
tive increase in the east cloud intensity. Future model 
calculations including the lo plasma torus are, however, 
definitely needed to clarify this situation. 

In summary, the model predicted changes in the shape 
and spatial orientation of the sodium cloud envelope and in 
the intensity distributions within the evelope, provide a 
simple explanation for the east-west shape and intensity 
asymmetries. These results follow directly from the lack 
of circular symmetry about Jupiter of the forces acting on 
the cloud atoms. The lack of this circular symmetry is 
introduced by the presence of the solar radaition force in 
the model calculations. The identification of the east-west 
intensity asymmetry in the model calculations with sodium 
atoms emitted within a somewhat narrow range of the veloc- 
ities near the satellite escape speed is significant. Care- 
ful model inversion of the east-west intensity data should 
then be sensitive in accurately determining relative weights 
in the initicl emission velocity distribution for the con- 
tributing velocity components in the 1.8 km sec ^ to 2.2 km 
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sec ^ velocity interval. This analysis would be somewhat 
complementary to the inversion of the two-dimensional inten- 
sity data for the east-west cloud shape asymmetry seen on 
the sky plane, which is primarily sensitive to the emission 
velocities between 2 km sec ^ and 3 km sec and to the 
line profile data of the sodium cloud, which is particularly 
sensitive to higher emission velocities, up to 15 km sec ^ 
or 20 km sec A consistent and simultaneous analysis of 
all three types of sodium data should provide the necessary 
overlapping constraints to deduce the real initial velocity 
distribution function for sodium atoms escaping from lo. 
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FIGURE CAPTIONS 


Figure 1 


Figure 2 


Figure 3 


East-West Changes in the Sodium Cloud. Calcu- 
lations performed in Paper I for the solar 
radiation perturbed cloud shape, indicated 
by the shaded area, and for the unperturbed 
cloud shape, indicated by the dashed boundary, 
are compared for diametrically opposite satel- 
lite phase angles. The shape of the cloud is 
the envelope formed by atoms, emitted from lo in 
the satellite plane, after 20 hours of flight 
time. 


Dimensionless Coordinate System for the Modified 
Circular Restricted Three-Body Problem in lo's 
Orbital Plane. The locations of the five classi- 
cal Lagrange points at L^, L 2 f L^ and L^, the 
satellite at x^, and the planet at X 2 in the (x,y)- 
plane are indicated. 

Sodium Electron Impact Ionization Lifetime in 
the lo Plasma Torus. The lifetime calculation 
is based upon the measured cross-section by 
McFarland (1965) , McFarland and Kinney (1965) 
and Zapesochnyi and Aleksakhin (1968) . For 
the plasma data, referenced in che text, the 
two-d-.mensional common temperature model of 



Bagenal and Sullivan (1981) v/as selected to 
describe the electron density. 

Figure 4 Interaction of the Sodium Cloud and the Plasma 
Torus. The spatially-projected overlap in the 
satellite plane of the lo plasma torus with the 
solar radiation perturbed cloud shape (solid 
boundary with heavy shading) and with the unper- 
turbed cloud shape (heavy dashed boundary) of 
Figure 1, is indicated for the diametrically 
opposite satellite phase angles. The cloud shape 
illustrates the envelope of sodium atoms after 20 
hours of flight time. 


Figure 5 Relative D 2 Intensity in the Bergstralh Slit- 

The relative contribution to ^.he solar resonance 
scattered D 2 intensity seen in the observing slit 
of Bergstralh et al. (1975, 1977) as determined 
by monoenergetic model calculations for the sodium 
cloud, excluding the effects of solar radiation 
pressure, is shown as a function of the emission 
velocity of sodium atoms from the inner exobase 
(2600 km radius) of lo for a satellite phase angle 
of 77.6 degrees- For each emission velocity, the 
same number of atoms were emitted from the exobase 
and a simple constant cut off lifetime of 20 hours 
was assumed. 
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Figure 6 Effects of Solar Radiation Pressure in the East 
Cloud. Calculated D 2 intensity contours of the 
sodium cloud are compared, with and without the 
effects of solar radiation pressure, for the 
indicated phase and sun angles. The outer dashed 
lines are for reference in comparing the left and 
right boundaries, and the center dashed line is 
for vertical .'lignment of the satellite location. 


Figure 7 Effects of Solar Radiation Pressure in the West 
Cloud. See legend to Figure 6. 

Figure 8 East-West Comparison of the Sodium Ci^oud. Model 
calculations for the D 2 intensity contours of 
the sodium cloud, excluding the effects of solar 
radiation pressure, are shown, for diametrically 
opposite satellite phase angles east and west of 
Jupiter. No east-west phase lag asymmetry is 
evident. See text for discussion. 

Figure 9 East-West Comparison of the Sodium Cloud. Model 
calculations for the D 2 intensity contours of 
the sodium cloud, including effects of solar 
radiation pressure, are shown for diametrically 
opposite satellite phase angles east and west of 
Jupiter. A distinct east-west phase lag asymmetry 
is evident. See text for discussion. 



Figure 10 Dependence of the Sodium Cloud on Emission 

Velocity. Model calculated D 2 intensity con- 
tours for the sodium cloud, including solar 
radiation pressure effects, are shown for 
diametrically opposite satellite phase angles 
and compared for emission velocities of 2.2 km 
sec 2.4 km sec ^ and 2.6 km sec See text 
for discussion. 

Figure 11 Dependence of the Sodium Cloud on Emission 
Velocity. See legend to Figure 10. 

Figure 12 Observational Slit for the Intensity Asymmetry 
Measurements. The rectangular observational 
slit of Bergstralh et al. (1975, 1977) , centered 
on lo, is shown in relation to the sodium cloud, 
the satellite orbit and the planet. 

Figure 13 Comparison of Observation and Model Results for 
the East-West Intensity Asymmetry. The obser- 
vations of Bergstralh et al. (1975) are compared 
to model calculations for the relative D 2 inten- 
sity seen through their rectangular viewing slit. 
Model calculations for four different monoener- 
getic values of the emission velocity are shown. 
See text for discussion. 
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Figure 14 Model Results for thti East-West Intensity Ratio. 


The east to west D 2 intensity ratio of Figure 13 
is shown as a function of the east to west phase 
angles for emission velocities of 2.0 km sec ^ 
and 2 . 1 km sec ^ - 

Figure 15 Lagrange Window in lo's Orbital Plane for Atoms 
Emitted from the Exobase of the Satellite. lo, 
surrounded by its exobase (the dashed circle) , is 
shown as depicted in an original illustration of 
Smyth and McElroy (1977). The radial coordinate 
indicates the flight time of an atom emitted radi- 
ally from the exosphere at the angular orientation 
about the satellite shown. Upon capture by lo, 
the 5° angular-resolution radial sector contain- 
ing that orbit is blackened for all later time. 
Ballistic orbits are indicated for a capture 
time of 50 hr. For greater capture times, black- 
ening indicates that atoms, initially on an escape 
trajectory, have been recaptured because of a close 
encounter with lo. 

Figure 16 Solar Radiation Pressure Modulation of the Sodium 
Atom Escaping from lo. The percentage of sodium 
atoms not escaping lo, which are emitted uniformly 
and at a constant rate from the inner hemisphere 
of the satellite exobas^ with monoenergetic emis- 
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Sion velocities of 2.0 km sec ^ and 2.1 km sec”^, 
are compared as a function of the lo phase angle 
for model calculation excluding (dashed curve) 
and including (solid curve) the effects of solar 
radiation pressure. See text for discussion. 
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